Quantum squeezing by a parametric resonance in a SQUID 
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We study rotating squeezed quantum states created by a parametric resonance in an open har- 
monic system. As a specific realization of the phenomenon we study a mesoscopic SQUID loop where 
the state preparation procedure is simple in principle and feasible with currently available experi- 
mental methods. By solving dynamics and calculating spectral properties we show that quantum 
fluctuations of SQUID observables can be reduced below their groundstate value. The measurement 
is introduced by coupling the SQUID to a transmission line carrying the radiation to a secondary 
measurement device. Besides the theoretical interest, our studies are motivated by an opportunity 
for a practical quantum noise engineering. 
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I. INTRODUCTION 

At the heart of the quantum theory lies the fundamen- 
tal principle of describing physically observable quanti- 
ties as Hermitian operators acting on Hilbert space of 
quantum states. Generally these operators do not com- 
mute, a trait giving rise to the fundamental uncertainty 
principle first formulated by Heisenberg. The uncertainty 
principle characterizes the statistical spread of the distri- 
butions corresponding to two different observables. If 
the operators representing the observables do not com- 
mute, the statistical variations of their observed distri- 
butions, frequently called uncertainties, cannot generally 
be arbitrarily small in a given state. However, the sta- 
tistical variation of a single observable is not limited in 
any way by the uncertainty principle. For two given 
non-commuting observables, a state in which the lower 
limit of the uncertainty principle is reached is called a 
minimum uncertainty state. Decreasing the uncertainty 
with respect to one observable results in the increase 
of the uncertainty of the other. The transfer of uncer- 
tainty from one observable to another with respect to the 
minimum uncertainty values is referred to as squeezing. 
The squeezing of quantum fluctuations was first stud- 
ied and experimentally verified in quantum optics, where 
the components of quantized electric field served as the 
squeezed observables Since then the phenomenon has 
been observed in superconducting circuits;^ and more 
recently, there has been a promising efforts to realize the 
squeezing in nanomechanical structures^ 

Due to the technological advances, the development in 
mesoscopic physics has been rapid in last decades, and 
many interesting quantum phenomena have been experi- 
mentally verified for the first time. Also a number of new 
and important measuring devices, necessary to observe 
elementary quantum phenomena, have been invented. 
In this paper we study the phenomena of the quantum 
squeezing in a mesoscopic Superconducting QUantum In- 
terference Device (SQUID), with emphasis on the cre- 
ation and consequences of the squeezing. However, the 
squeezing is an example of a manipulation of quantum 
fluctuations which could have direct applications in a 



practical quantum measurement. 

Experimentally the squeezing of the quantum fluctu- 
ations has been demonstrated in microwave frequencies 
by constructing a Josephson parametric amplifler,'* with 
a 40% reduction of the vacuum noise. Theoretically the 
squeezing has been studied in a similar SQUID ring we 
consider hereA The squeezing mechanism we study is dif- 
ferent from one studied in Ref. 6 which was based on a 
rapid decrease of an external magnetic flux to switch on 
the Josephson coupling. We consider a parametric in- 
stability in a harmonic regime,'' a well-known procedure 
to create the squeezing and applicable in various differ- 
ent systemsi^ It can be realized with an elementary flux 
control in an rf SQUID. In the limit of a negligible dis- 
sipation the magnitude of the squeezing is exponential 
in short times and, it rotates between the charge and 
the magnetic flux of the SQUID. A strong dissipation 
compensates the resonant driving and leads to a rotating 
quasistationary state where uncertainties periodically go 
below their groundstate values. We calculate the noise 
spectrum of periodic squeezed states and discuss the mea- 
surement problem by analyzing a setup where the SQUID 
is coupled to a transmission line. 



II. SQUEEZING BY PARAMETRIC 
RESONANCE 



The parametric resonance is a well-known physical 
phenomena in a classical harmonic oscillator (HA)i^ A 
periodic time-dependent perturbation can lead to large 
effects, which are most signiflcant when the period of the 
perturbation is twice the resonance frequency. Driven in 
resonance, the energy of the system grows exponentially 
in time. This is true for a quantum HA also.^ In addition, 
starting the resonant driving from the groundstate of HA 
it results in a rapid squeezing of the uncertainties of the 
conjugate observables of HA. Consider the Hamiltonian 
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where p and x are the canonical variables [x,p] = ih, m 
is the mass constant and the time-dependent dependent 



perturbation is included in lo (t) 



' Acos2cuot. Fur- 



thermore, let us assume that the time-dependent part 
can be treated as a perturbation, \A\ <C ujq. The ground- 
state of the free Hamiltonian {A = 0) in the posi- 
tion representation is ipoix) = TT^^^^a^/^exp{—ax^/2), 
where a = 1//q = mujo/h is an inverse square of the 
characteristic zero point length scale. An approximate 
solution to the Schrodinger equation corresponding to 
the Hamiltonian of Eq. |^ and the initial condition 
ip{x,t = 0) = ipoix) can be obtained analytically in the 
Gaussian form 

^x, t) = {^^\ exp {-ait) x^l2) , (2) 



where 



a{t) 



i (1 — ie^)sin(tJo^) + (1 + ie «)cos(woO 
l'^ -(1 + ie2?)sin(wo<) + (1 - ie2«)cos(wot) 



(3) 



with ^ = At/Aojo- Physically Eq. |5J describes a wave 
packet centered at x = 0, but whose width oscillates in 
time with an exponentially growing amplitude. When the 
statistical spread in the real space is the widest, the cor- 
responding spread in the momentum space is the small- 
est. Defining the dimensionless quantities x' = x/Iq and 
p' = plo/h we discover that in the state (O we have 
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|([a;',p'])|/2 = 1/2 and the 



state remains minimum uncertainty state at all times. 
The individual deviations may be written as 
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where /3= {Rea{t))-^/'^ /Iq. From Eq. Q we obtain 
Rea{t) — (sinh^^+ cosh^^ 

— 2 sin(2wot)sinh^cosh^ (5) 

The maximum uncertainty oscillates with a frequency 
2u!o between the position and momentum. At the times 
corresponding to sin(2ajoO = il the squeezing is given 
by /3 = e=F«. The solution 101 is very accurate as long 
as the condition A ujq holds, li A < O.Iojq the ana- 
lytical solution agrees excellently with the exact solution 
up to the squeezing factor 5, see Fig. ^ After that the 
lower value does not decrease further as in the analytical 
solution even though the general agreement is accurate. 

The above example illustrates the mechanism behind 
the squeezing scheme we will now apply to a flux- 
controlled rf SQUID loop described in Fig. [21 The Hamil- 
tonian for the system can be written as 
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E.,cos{^{t)2e/h)cos{(l)2e/h), (6) 



where C is the capacitance of the junctions, is the self- 
inductance of the loop and <I>(t) is the flux bias externally 



applied through the smaller loop.^^ The quantities Q and 
(j) are the canonical variables corresponding to the charge 
in the capacitor and the magnetic flux in the loop and 
satisfy the commutation relation [(f), Q] — ih. Thus the 
magnetic flux plays the role of the position coordinate 
and the charge Q corresponds to the momentum in the 
standard one-particle quantum mechanics. The flux value 
corresponding to the oscillator length Iq of the harmonic 
part of © is (/>o = Vj 



where = ^Ls/C. 




FIG. 1: (Color online) Uncertainty of the magnetic flux 
\/^A4>' in time. The parametric driving with a strength 
B — O.ltJo makes the uncertainty to oscillate with the fre- 
quency 2u!o below and above the groundstate value marked 
by the black dashed line. The blue line corresponds to the 
analytical solution, the red line to the exact numerical solu- 
tion and the green line to the exact solution with the AC flux 
modulation 

We treat the Josephson term as a perturbation and 
therefore consider the parameter range Ls <C Lj, where 
we have defined the Josephson inductance as Lj = 
h^/Ae^Ej. If the quadratic potential is strong, the flux 
particle is restricted to the linear regime of the SQUID 
and the term cos{(j)2e/h) may be expanded to the sec- 
ond order. Supposing that the condition for the validity 
of the approximation (jy^ < $o = h/2e is satisfled, the 
Hamiltonian 10 takes the form 
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cos($(t) 2elh) = 



where 



1 - ^cos($(t) 2elK) 



(7) 



(8) 



with luq — {LsC)~ . In the cosine expansion the small 
term constant in cj) has been dropped from Eq. ((T)). Now 
we have established a connection between Eqs. and 
Q provided that the external magnetic flux is modu- 
lated according to ^{t) = hujot/e. This could be achieved 
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by increasing the magnetic field linearly inside the con- 
trol loop. The temporal evolution of the SQUID start- 
ing from the groundstate is readily solved by identifying 
A ujq{Ls/Lj) and Iq 4>o and applying Eqs. Q 
and lO- The dimensionless operators (j>' — 4>/4'Q and 
Q' = Q(j)Q/h satisfy the minimum uncertainty condition 
A^'AQ' — 1/2 and the individual deviations behave as 
those in Eq. 

The linearly increasing magnetic flux is inconvenient 
practically since the magnetic field should be restricted 
to the control loop and it is difficult to isolate at high 
field values. In addition, the AC fields are obtained and 
manipulated more easily. From the approximate equiva- 
lence cos (2 x) « cos (3 sin(a;)) can be inferred that if the 
magnetic flux is modulated with an AC field according 
to 

= (3/27r)$osin(cc;oO, (9) 

then cos{^{t) 2e/h) « cos{2LUQt). Hence the magnetic 
fiux should be modulated at the resonance frequency such 
a way that its amplitude is roughly a half of the flux 
quantum. 

For later purposes it is convenient to introduce the 
second-quantized creation and annihilation operators in 
the photon number Fock space. The SQUID observables 
may be written as 

+ Q^^^y-a). (10) 

Operators a, satisfy the usual bosonic commutation 
relations. In this language the Hamiltonian Eq. Q is 
transformed to 

H = nwQ{a)a + i) + S cos (2 u}Qt){a + a) f, (11) 

where B = HujqL/ALj. In the case where expectation 
values of and Q vanish, as above, the root mean square 
deviations may be written as 

V 2 

AO' = ^ (-(a^) _ (at2) + (aat) + (a^a))'^^ . (12) 
V2 

III. COUPLING THE SQUID TO A 
TRANSMISSION LINE 

In this section we introduce a measuring scheme in 
which the SQUID is coupled to a transmission line (TL) 
which serves as a wave guide carrying the radiation away 
from the system. It also plays a role of a generic mea- 
surement device which causes a Markovian back-action 
to the measured system. In addition, the TL provides a 
practical theoretical model for studying the environmen- 
tal effects to the squeezing. 




FIG. 2: Schematic representtion of the SQUID loop. The 
crosses represent Josephson junctions and the physical quan- 
tities (j} and Q correspond to the magnetic flux through the 
main loop and the charge at the junctions. The flux $(t) 
through the smaller loop is controlled by an external mag- 
netic fleld. 

The study of a non-classical radiation requires a 
quantum-mechanical treatment of the SQUID and the 
TL. Quantum mechanics of a TL is essentially similar 
with a one-dimensional quantized electromagnetic field 
and it has been discussed in Ref. JTj Consider the sys- 




FIG. 3; Transmission line inductively coupled to the SQUID. 

tem in Fig. El where the SQUID is magnetically coupled 
to a TL. The magnetic flux of the nearby SQUID pene- 
trates into the TL. The inductance and capacitance per 
unit length in the TL are I and c and the wave velocity 
is given by w = 1/ Vic. 

The Lagrangian for the line can be written as 

where j{x, t) and q{x, t) are the local current and charge 
densities. Following Ref. 11 we introduce a new variable 



e{x,t) = / dx'q{x',t), (14) 
Jo 

which, together with the continuity equation, allows one 
to write the Lagrangian 
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From the total charge neutrality that it follows that 
0{O,t) = 0{L,t) = 0, so 9 can be expanded by the eigen- 
modes 



(16) 



where bk — kn/L and the Lagrangian takes the form 



(17) 



The quantization of the quadratic Lagrangian ()17|l is well 
known. The Hamiltonian can be written as 



H 



TL 



(18) 



where [^^,11^] — ih. This can be further written with 
the help of bosonic creation and annihilation operators 
Cfe and ct as 



fkukc 1 



{cu(t)+cl{t)) 



and the Hamiltonian (|18|l becomes 



1/2). 



(19) 



(20) 



which can be solved as 

Ck{t) =Cfc(0)e-^"^* + 

l^y^e— ^*^V(i')e^"''*'di', (23) 

where (j){t) is the flux operator of the SQUID. From this 
expression one can see that the TL problem can be solved 
as soon as the the solution to the associated SQUID prob- 
lem is known. The flux (l>{t) is the only SQUID variable 
entering to the expression and can be solved inde- 
pendently provided that the TL variables can be success- 
fully eliminated from the SQUID dynamics. Substituting 
the expression H23|l into the formula for the voltage of the 
TL 



Vix,t)=Y^^^cos'^[c,(t)+ci{t)], 

k 



we get 



V{x, t) = Vo{x, t) + ^^^{t - x/v), 



(24) 



(25) 



where Vo{x, t) is the voltage operator of the free TL and 
the second term is proportional to the retarded SQUID 
field. Deriving Eq. ()25|l we have assumed that the level 
separation in the TL is much smaller than wq, so the 
TL will act as a waveguide and not as a resonator. The 
current of the TL has a similar expression. 



The eigenfrequencies are defined as LOk = PkV- The volt- 
age and current operators in the TL can be derived from 
9 as Vix,t) = {l/c)da:9{x,t) and I{x,t) = 9{x,t). The 
magnetic interaction between the SQUID and the TL 
couples the operators Is — 4'/Ls and I{x,t), where Is 
and Ls are the current and the inductance of the SQUID 
loop. Supposing that the SQUID couples to each mode 
approximately as strongly, the interaction may be writ- 
ten as 



hujk , 



4)' 



(21) 



where M is the mutual inductance between the SQUID 
and the field mode. 

Now we could proceed to two directions. On one hand 
we could study the SQUID Hamiltonian by trying to 
eliminate the transmission line variables (and other possi- 
ble environment modes) . On the other hand we are inter- 
ested in the TL output radiation to measure the SQUID. 
The first problem leads to a master equation description 
of the SQUID and it will be studied in the next section. 
The latter problem will be pursued here. 

From the Heisenberg equation of motion for the oper- 
ators Cfc we obtain 



OtCk = -iUJkCk + -;—\ TTl 



hLV 



(22) 



IV. THE ROLE OF THE DISSIPATION AND 
DECOHERENCE 

A. Dissipative state evolution of the SQUID 

So far we have treated the SQUID as an ideal system 
without any dissipation and decoherence effects. Real 
mesoscopic systems are much larger than the systems 
usually studied in atomic and molecular physics and gen- 
erally need to be considered as open quantum systems. 
Also an important source of decoherence is the quantum 
measurement process in which the system is necessarily 
coupled to the external world. A state of an open quan- 
tum system can no longer be represented by simply a 
state vector making it necessary to introduce the density 
operator. The dissipation, energy absorption and sup- 
pression of quantum coherence effects can be naturally 
discussed within the density operator approach. 

Assuming that the level separation in the SQUID loop 
is significantly higher than the temperature of the en- 
vironment hujQ > fceT, the dominating environmental 
effect is the spontaneous emission. The dynamics of the 
density operator after the Born-Markov elimination of 
the electromagnetic field modes can be described by a 
Lindblad master equationii^ of the form 

dtp = — — [H, p] + K{2apa) — a^ap — pa'a), (26) 
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where the first term is responsible for the Hamiltonian 
evolution by the operator Ullfl and the second term de- 
scribes the spontaneous emission into the TL, which is 
assumed not to bring any signal back to the SQUID. 
The coefBcicnt k is related to the quality factor of the 
circuit by K = o^o/Q- In the case of a free SQUID the 
Q- factor can be as high as 10^ — 10^. The coupling to 
the TL should decrease Q significantly to guarantee an 
efficient measurement of the SQUID. The value of kappa 
can be estimated then duS k [M / LsY Zq/ Z'^i^, where 
^TL — \/W^- The disadvantage of decreased Q is the 
increased decoherence which is hostile to the squeezing 
process. 

B. Expectation values 

Starting from Eq. H26|l it is possible to derive the equa- 
tions of motion for the expectation values of the creation 
and annihilation operators, 

9t(a) — {—iuJo — K — 2iBcos ( 2uJot)) (a) 

- 2iBcos{2LUot){a^), 
9t(a^) = {+iuja — K + 2iBcos ( 2(^0^)) 

+ 2iBcos{2LUot){a), (27) 

as well as for those of the quadratic operators 

dt{a^a) = 2iBcos{2uJot) ((a^) - (a^^)) - 2^(0^0) 
dt{a'^) = — 2icjo(a^) — iBcos{2ijjQt)'K 

(4(a2) -I- 4(ata) + 2) -2Hi{a'^) 
aj(at2) ^ 2iwo(a^^> +iBcos(2uJot)x 

(4(at2) + 4(ata) + 2) - 2n{a^'') . (28) 

These two present a closed set of equations. The missing 
combination (aa^) can be deduced from the commutation 
relation [a, a^] = 1. By solving the above equation system 
we obtain the uncertainties (|12|l . The nature of solutions 
to the system H28(l depends on the relative magnitudes 
of K and B. If k > B, the dissipation will eventually 
compensate the resonant driving and the solutions are 
periodic and bounded. Otherwise the expectation values 
grow increasingly in time. This, however, does not imply 
increasing squeezing of quantities (|12|l . In the presence 
of dissipation the lowest value of the fluctuations has a 
finite lower limit. 

By expanding the expectation values in Eq. H28(l to a 
Fourier series and truncating them to the fourth order, 
we found an accurate analytical form for the periodic 
solution. In the limit k — s- i? -I- 0, the lower limit of the 
periodic squeezing is about 0.75 times the vacuum value 
of A(j)' and AQ', see Fig. El 

The nature of squeezed states in the presence of the 
dissipation is somewhat different from the pure squeezed 
states. The pure squeezed states are linear combina- 
tions of the even photon number states. The dissipation 
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FIG. 4: (Color online) Two periods of the stationary solution 
K = 1.5B. The blue curve is Atf)' and the red curve {a' a). 
The yellow and the green curves correspond to Re{a^) and 
Im (a^) respectively. The black horizontal dashed line marks 
the ground state value of Ad)' . 
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FIG. 5: (Color online) Uncertainty Ac/)' of the bounded peri- 
odic solutions K = 1.5B (blue), k = 1.3B (red) and k = 1.1_B 
(green). The black horizontal dashed line marks the ground- 
state value of A(/!>'. The minimum of the squeezing in the 
periodic solutions approaches to lower bound of about 0.75 
times the groundstate value. 

changes the spectral content so that the odd number pho- 
ton states are also occupied. The resulting state is no 
longer a minimum certainty state, see Fig. El 

C. Two-time correlation functions 

In order to find the spectral properties of the signal 
radiated into the TL, we also need to evaluate two-time 
correlation functions, such as (0^(^2)0(^1))- According 
to the quantum regression formulaic, the function pair 
{A{t)a{t')) and {A{t)a'' {t')) obey the same differential 
equations than {a{t')) and (a^(i')) for arbitrary opera- 
tor A{t). Choosing A{t) to a{t) we get 

Of (a(t)a(t')) = {-iiuo - k - 2iBcos ( 2cjo<')) 
- 2iBcos ( 2LUot'){a{t)a''{t')), 
df {a(t)a'^ {t')) = {+iLUo - k + 2iBcos ( 2ujot')) {a'' {t)a{t')) 
+ 2iBcos ( 2LUot'){a{t)a{t')) (29) 

for t' > t. A similar pair of equations hold for {a'' {t)a{t')) 
and (at(t)a^(t')). Solving Eqs. (|29|) and the correspond- 






FIG. 6: (Color online) Evolution of the uncertainty product 
A</)'A(5' from the groundstate to the periodic steady state 
K = 1.5B. After reaching the steady state, the uncertainty 
product oscillates with a frequency 4ujo about the constant 
value. 



ing equations to the other pair, we can deduce an arbi- 
trary second-order correlator in the whole time domain 
by using symmetries of the correlators. 
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FIG. 7: (Color online) Wigner function of the groundstate 
(left) and the periodic squeezed state k = 1.55 (right) in 
the (0', (5')-plane. The squeezed state rotates clockwise as 
indicated by the arrows. The ellipse makes 2tt rotation in 
time T — 27r/cJo. 



A. Time-averaged spectrum 



D. Wigner function of the squeezed state 

An intuitive way to perceive the squeezing phenomena 
is to plot the Wigner function 

pw{x',p') = (2^)-i J {x'~-y\p\x'+-y)e'P ^ dy. (30) 

The Wigner function is one of the quantum-mechanical 
analogues to the phase space probability distribution. 
Even though it is not a genuine 2-dimensional prob- 
ability distribution, it has the property of giving the 
one-dimensional marginal distributions in arbitrary di- 
rections in the phase space {x',p') and can therefore be 
used directly to illustrate the uncertainties of observables. 
The groundstate of the SQUID has a circular symmetry 
which is distorted to an ellipse by squeezing, see Fig. 
[7| For ideal squeezed states the principal axes of ellipti- 
cal contour lines are inversely proportional reflecting the 
minimum uncertainty property. The primary effect of 
the dissipation is to broaden the shorter axis resulting in 
an excess uncertainty. 



V. SPECTRUM OF PARAMETRIC 
RESONANCE AND SQUEEZING 

Based on Eq. the spectral properties of V{X, t) 

are directly related to spectral properties of 4> and the 
free TL spectrum. Here we evaluate (/)-noise and, at the 
end of the section, discuss parameters relevant to the TL 
output noise. 



The spectrum contains information about radiation 
emitted from the studied system. In Ref. ^3 an excitation 
of a two-level system with an energy splitting = AE 
is considered as a spectrum analyzer. The probability 
to detect the analyzer system in either of its states is 
proportional to 

f f < A{ti)A{t2) > e-'-(*^-*i)dtidi2 (31) 
Jo Jo 

at io = ±r2, where A is an operator of the studied noise 
source which couples to the analyzer two-level system. 
Supposing that the system is in a steady state and the 
correlation time is short, the expression (|31|l reduces to 
i5'A(w), where 

/oo 
< A{t)A{0) > e"^^dT. (32) 
-oo 

The expression (|32f) is positive definite and determines 
the transition rates of the analyzer system. In the strong 
damping regime the SQUID radiation is periodic and the 
steady-state expression H32|l is not appropriate as such. 
Rather, the detector rates are proportional to the time- 
averaged noise power 

Sa{uj)= f ^[ <A{T + t')A{t')>e'^^dt'dT, {33) 

J -oo ^ Jo 

where T is the period of the motion. The expression H33|l 
generalizes (|32|) in the sense that so defined spectrum is 
positive definite and gives excitation probabilities of the 
analyzer system in the long-time limit t ^ T — I-kIloq. 
Since the SQUID couples to the TL with the operator 
= (a -I- a^^)/\/2, the (/)'-spectrum is relevant in the TL 
radiation. 
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From the time-averaged observables it is difficult to 
detect the squeezing directly. This can be understood 
by considering the Wigner function. The fast rotation 
of the elliptic figure in Fig. [7| averages out to a circular 
shape where no squeezing is present. The (/)-spectrum 
of the SQUID in parametric resonance resembles to 
a spectrum of a dissipative harmonic oscillator at finite 
temperature, see Fig. |S1 The spectrum displays peaks at 
frequencies —luq and loq corresponding to emission and 
absorption processes. A study of a fine structure also 
reveals weak resonance peaks at frequencies —Sloq and 
SajQ. These are the higher harmonics produced by the 
parametric driving process. In principle there is also reso- 
nance at higher odd multiples of co but they are extremely 
weak compared to —loq and loq peaks. The distance be- 
tween the peaks is 2a; which characterizes the periodicity 
of the rotating squeezing. 
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FIG. 8: (Color online) (a) Time-averaged noise spectrum 
S^i{ui) of the periodic squeezed state k = O.lSoJo = 1.255. 
The Spectrum has resonances at frequencies ±ti;o. (b) The 
fine structure of spectrum (a) reveals the weak resonance 
peaks at ±3ljo. 



rotation of the harmonic creation and annihilation oper- 
ators. In the new operators the orientation of the squeez- 
ing is almost static and determined by the angle 0. The 
Wigner functions in Fig. appear frozen in time. Only 
a tiny time-dependent deformations resulting from the 
driving remains. The noise varies periodically as a func- 
tion of Q as the orientation of the squeezing is rotated, 
see Fig. O For a finite range of d the SQUID noise goes 
below its groundstate value which is a direct evidence of 
reduction of the quantum fiuctuations. For the parame- 
ters of Fig. El the groundstate noise at zero frequency is 
approximately 2.5 times higher than the minimum noise 
corresponding to = 0. The noise curves of squeezed 
states are slightly distorted respect to the origin, which 
indicates that there is a tiny asymmetry respect to the 
half axes of the ellipse in Fig. A parametric driv- 
ing in the presence of a dissipation does not lead to a 
completely symmetric squeezing. 

A phase-sensitive measurements can be realized, for 
example, with a help of a parametric amplification 
processr'^'^ A harmonic modulation of the measurement 
signal with the frequency can be used to measure the 
rotating spectrumii^ 
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FIG. 9: (Color online) Noise spectrum S^i (w) of the squeezed 
state K = 0.15ti;o = 1.55 in the rotating operators. The 
resonance peaks are now moved to the origin. The initial time 
is chosen so that the minimum noise curve (blue) corresponds 
to ^ = 0. The green, red and yellow curves correspond to the 
orientations d = 7r/4, 7r/2, 37r/4 respectively. The black curve 
represents the groundstate noise. The minimum noise curve 
remain below the groundstate noise, which is an evidence of 
squeezing. 



B. Spectrum of rotating operators 

The detailed analysis of squeezing can be performed in 
a phase-sensitive measurement. For that purpose we de- 
fine the rotating frame operators as b(t) — exp{iLUQt + 
e)a{t) and h^{t) = exp{-iuJot - e)a^{t). The rotat- 
ing SQUID observables are defined analogously as 0' — 
{b{t) + b^t))/V2 and Q' = i{b^ (t) - b{t))/V2. The rotat- 
ing phase factors in b{t) and 6^(t) compensate the natural 



C. Transmission line output 

The TL observables are sum of a free line observables 
and a retarded SQUID radiation contribution, see Eq. 
(|25|l . The voltage noise of the TL is of the form ^^(a;) = 
Sy + g'^Sfj,! {ijj), where Sy is the vacuum noise of the TL. 
The coupling constant g depends on the actual material 
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parameters and can be deduced from Eq. (|25|l : 

where = Lg/C. For the phenomenon to be mea- 
surable, the vacuum noise of the TL should not over- 
whelm the SQUID noise. From Eq. |Mlwe can estimate 
the order of magnitude as g'^Scj,>{LOo) ~ {M / LsY'hujQZ^. 
At low temperatures the free TL vacuum noise is 
Sy = HlooZti^, where Ztl = \/lJc is the characteris- 
tic impedance of the TL. Then the relative effect is of 
order g^S^>{u;o)/S^{uJa) = {M / Lsf Zo/ Z^i^. For realis- 
tic values M/Ls = 0.2, ^ 50017 and Ztl = 5017 the 
fraction becomes 5^5*0' (t^o)/ ■5'$/ (^o) = 0.4, which sug- 
gests that the SQUID noise is significant contribution to 
the TL output radiation in favorable conditions. 

VI. CONCLUSIONS 

The parametric harmonic driving creates rotating 
squeezed quantum states in the SQUID ring. If the 



environmental damping of the SQUID is negligible the 
squeezing of uncertainties is exponential for short times. 
In the presence of a strong damping the magnitude of 
the squeezing is stationary. The minimum uncertainty 
in the flux and charge go below the groundstatc value 
periodically. The phenomenon enables a quantum noise 
engineering which plays an increasingly important role in 
the quantum measurement theory as well as in the de- 
sign of practical quantum devices. In principle, squeezed 
harmonic systems could be used as an ultra sensitive 
measurement devices of external perturbations coupled 
to it.^ 



The experimental creation of squeezed quantum states 
in a SQUID by a parametric driving is feasible with the 
current experimental methods. We calculated the rele- 
vant noise properties of the periodic squeezed state. By 
introducing a coupling to the transmission line, we an- 
alyzed the fully quantum mechanical emission spectrum 
of the SQUID and discussed briefly the conditions of an 
experimental verification of the phenomenon. 
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